| |
INTRODUCCIÓN
La resolución de los picos obtenidos en técnicas separativas, como la electroforesis y la cromatografía, depende de muchos factores como la longitud y calidad de las columnas utilizadas, las fases móviles y estacionarias de las mismas, los electrolitos de fondo y tensiones de polarización usados en el caso de la electroforesis y muchos otros que sería prolijo enumerar.
Siempre se busca elegir las variables instrumentales que optimizan la resolución y los tiempos de elución para los analitos de mayor interés, pero es inevitable que algunas veces aparezcan en los registros otros pìcos que se superponen o solapan con alguno de los picos que interesa cuantificar, lo que obligaría a modificar las condiciones de trabajo para evitar ese solapamiento o bien disponer de algún algoritmo que permita descomponer el pico obtenido en sus componentes, para asi ser capaces de cuantificar correctamente los analitos de interés que formen parte del pico compuesto.
Esta última opción es la elección correcta cuando no se quieren modificar las condiciones instrumentales que se han elegido y que funcionan adecuadamente para la mayoría de los análisis que se realizan.
La dificultad para aplicar un algoritmo de este tipo dependerá de varios factores, pero entre ellos y de mucha importancia estará la proximidad entre los tiempos de elución de los picos que no llegan a resolverse, como puede verse en la Figura-0 donde se presenta la simulación de tres electroferogramas con tres picos que eluyen a dos tiempos fijos el primero y segundo y el tercer pico con una diferencia de tiempo de 20s, 15s y 10s, respectivamente, con respecto al segundo pico.
Se ve en la figura, que en el primer y segundo electroferogramas la cuantificación se puede realizar con cierta facilidad, mientras que en el tercero, aun siendo posible en este caso, es mas difícil dada la proximidad de los tiempos de elución de los picos segundo y tercero.
Si los tiempos de elución de los dos picos coincidieran exactamente no se tendría opción para diferenciarlos ya que solo se observaría un único pico resultado de la composición de los dos, excepto en el caso de que las características de ambos picos fueran significativamente diferentes entre sí y habría que recurrir a modificar alguna de las condiciones de trabajo que permitiera una separación mayor de los tiempos de elución de los dos picos.
En tal caso solo la experiencia del analista y las características del análisis le permitirán detectar que existen dos picos superpuestos en el mismo tiempo de elución.
Algunos de los algoritmos desarrollados para la descomposición de picos compuestos utilizan diferentes tipos de funciones matemáticas para representar los picos dependiendo de las características de estos, funciones gaussianas y otras que se ajusten a la forma de los picos para obtener por mínimos cuadrados un ajuste de la suma de varias de estas funciones a la curva del pico obtenido experimentalmente. Alguno de estos algoritmos está disponible en Internet (1).
Son muchas las funciones disponibles que se han ido incorporando en las nuevas versiones, dada la dificultad de representar las diferentes variantes de picos experimentales que se obtienen en los análisis.
Una de tales dificultades, que se presenta con cierta frecuencia, es que los picos obtenidos
experimentalmente no son simétricos con respecto al máximo o mínimo del pico, sino que presentan diferentes tiempos de subida y bajada y aunque existe una serie de recomendaciones (2,3) para evitar esta asimetría no siempre es posible evitarla.
Por ello en este texto se intenta una aproximación diferente, con un solo tipo de función matemática para representar los picos pero con unas pocas variables que permiten una gran diversidad en la forma de los picos simulados, cuya suma habrá que ajustar al pico experimental.
Como punto de partida se considera que los picos no tienen por que ser simétricos con respecto al máximo o mínimo, aunque también se puede con las variables adecuadas conseguir que lo sean si es necesario.
Este aspecto se presenta algo exagerado en la Figura-1 en la que se presenta un pico fuertemente asimétrico con respecto al tiempo en que se produce el mínimo de intensidad del pico, que se denomina t0 y cuya intensidad en el mínimo se denomina Cm.
La linea de base se considera plana y constante, con valor de intensidad Cb y los tiempos de bajada y de subida se representan por Tb y Ts, respectivamente.
Con estas variables se construye la función matemática que representa el pico, tal y como se describe a continuación. |  |
Un pico de esa forma se puede representar matemáticamente por la expresión : C( t ) = Cb . e-g( t )
donde Cb es la intensidad en la base del pico y g( t ) es tal que : g( t ) = 0 para t 0 + Ts < t < t0 - T b
y g( t ) = Km . [ Tb + F . ( t - t0 ) ] para t 0 - Tb <= t <= t0 + Ts , donde F=1 para t 0 - Tb <= t <= t0 y
F = - Tb / Ts para t0 < t <= t0 + Ts , siendo F = -1 cuando el pico es simétrico con respecto a t0 .
En la expresión de g ( t ) el parámetro Km afecta a la intensidad del mínimo del pico, ya que :
para t = t0 es g ( t0 ) = Km . Tb y C( t0 ) = Cb . e -Km . Tb y como a C( t0 ) se le ha llamado Cm ,
se tendrá que : Cb / Cm = e Km . Tb , de donde Km = ( 1 / Tb ) . ln ( Cb / Cm ) .
De este modo, a partir de esta expresión para representar la forma del pico se ve que los parámetros para caracterizar el pico son ; a) El tiempo de elución en el mínimo ( t0 ) . b) El tiempo de bajada hasta el mínimo ( Tb ) . c) La constante de intensidad ( Km ) . d) La constante ( F ) que relaciona el tiempo de bajada con el de subida, y e) La intensidad en la base ( Cb ).
REFERENCIAS
1.- blog.originlab.com/how-do-i-perform-peak-deconvolution/
2.- secyta.es/es/es/node/25
3.- waters.com/nextgen/es/es/education/primers/troubleshooting-guide-for-hplc.html
AJUSTE DE DOS PICOS
El caso mas frecuente que suele presentarse en la práctica es la superposición de dos picos en el registro del cromatograma o electroferograma y dependiendo del grado de resolución de los picos, resultará mas o menos sencilla la descomposición del pico compuesto para obtener y poder cuantificar los dos picos individuales incluidos dentro del pico obtenido en el registro.
Para optimizar el ajuste de los dos picos al pico experimental se traslada éste con la escala de tiempos expandida, con la que también se generan los picos que al combinarse obtendrán la mejor aproximación al pico experimental y en la que pueden ser cuantificados, con referencia a patrones previamente registrados.
Como ejemplo para este caso de ajuste de dos picos se toma un pico que se obtiene en electroforesis capilar en placa con detección por conductividad, usando como electrolito de fondo ácido acetilsalicílico, una columna de 10cm de longitud, una inyección de 2ul y una tensión de polarización de 500v.
Con estas condiciones se obtiene una buena separación de K, Na, Li y Ca, pero la presencia de Mg interfiere con el pico del Na debido a que tienen tiempos de elución muy similares, separados por aproximadamente 10s y esto imposibilita la cuantificación directa de cualquiera de estos dos metales a partir del pico obtenido experimentalmente.
Para ello se ha realizado la electroforesis de un jarabe que contiene solo Mg y Na y también se han inyectado patrones de K y Na. El fabricante asegura que el jarabe contiene 100 mg/ml de Mg y que podría contener hasta 23 mg/ml de Na.
Dada la viscosidad del jarabe se diluye en medio ácido a 1:10000 para extraer los metales presentes y facilitar la inyección y el registro dentro de escala de los picos obtenidos.
En el registro obtenido se observa que efectivamente no hay presencia de K ni Ca ni Li y solo aparece un pico, compuesto de Mg y Na presumiblemente, en el que interesaría cuantificar la concentración de Na, usando el patrón de Na inyectado, para ver si cumple con lo especificado en el prospecto del producto.
Se traslada este pico, obtenido experimentalmente, con una escala de tiempo expandida como puede verse en la Figura-2, en la que se observa la fuerte asimetría del pico y en cuya escala se van a construir los picos que han de combinarse para igualar al obtenido en el registro experimental.
El pico registrado tiene una base plana y de valor constante por lo que se ajusta al valor Cb para la construcción de los picos. |
Figura - 2 |
Una vez trasladado el pico experimental a la escala de tiempos, donde la abscisa total es ahora 10 min es decir expandida 6 veces con respecto a lo indicado en la leyenda, hay que tener en cuenta que los tiempos no son los que se leen en el eje de abscisas de la figura, sino tiempos relativos al mínimo del pico con un total de 360 s entre comienzo y final del registro trasladado, por lo que la construcción de los picos tendrá que hacerse con arreglo a estos límites de tiempo.
El primer pico que se construye matemáticamente tendrá que cumplir con los requisitos de tiempo de bajada, tiempo del mínimo e intensidad, los cuales se fijan inicialmente para ajustarse a los del pico experimental .
Se puede ver en la Figura-3 que las curvas de bajada se ajustan adecuadamente y el valor de Km podrá quizá ser modificado posteriormente, en caso de que sea necesario. |
Figura - 3 |
No se conoce inicialmente cual puede ser la diferencia de tiempos entre los mínimos de los dos picos superpuestos, pero se puede tener una estimación a partir del conocimiento de los tiempos individuales de elución, en caso de que se conozcan los componentes como sucede aquí.
Se sabe que esta diferencia de tiempo suele estar en el entorno de 10 s, por lo que se podrá ensayar dicho tiempo para el mínimo del segundo pico y lo que es una desventaja, que es la proximidad de los tiempos de elución, se convierte ahora en una ventaja ya que los tiempos de bajada dependen de los tiempos de elución, que en este caso al ser tan próximos permiten suponer con cierta garantía que dichos tiempos de bajada serán iguales.
Al asignar el mismo tiempo de bajada y un tiempo al mínimo desplazado en 10 s con respecto al primer pico y asignando un cierto valor al Km del segundo pico, se observa que la curva obtenida se aproxima mas a la experimental.
Se puede ver en la Figura-4 que aún queda por asignar un valor a la relación entre tiempos de bajada y de subida que podrá mejorar el ajuste de las curvas. |
Figura - 4 |
En el caso de no conocer los componentes del pico y por tanto no poder hacer una estimación de la diferencia de los tiempos de elución, habría que hacer un proceso de prueba y error para conseguir el mejor ajuste posible.
Dada la proximidad de los tiempos de elución y las características de los picos de ambos componentes, en este caso se puede considerar que la constante que relaciona el tiempo de bajada con el de subida es la misma para ambos picos.
Asignando a F un valor que haga coincidir las pendientes de las curvas de subida en los picos experimental y el simulado, se observa que la curva obtenida para el pico simulado se aproxima mas a la del experimental.
Se puede ver en la Figura-5 que aún puede mejorarse ligeramente el ajuste de las curvas, retocando algunas de las constantes asignadas, en concreto los tiempos de elución y las constantes de intensidad de los picos. |
Figura - 5 |
En este momento, pequeñas variaciones de alguna de las constantes elegidas podrán producir variaciones significativas en el ajuste de las curvas de los picos experimental y simulado, que mejorarán o empeorarán dicho ajuste pero que solo producirán efectos del orden de aproximadamente un 10% en los valores resultantes de la cuantificación de los picos, no obstante lo cual es necesario intentar el mejor ajuste posible.
En concreto, desplazando 2s el tiempo de elución del segundo pico hasta los 12s con respecto al tiempo del primero y retocando un poco las constantes de intensidad para ambos picos se consigue mejorar ligeramente el ajuste del pico simulado al experimental, como puede verse en la Figura-6.
En ella se incluye la cuantificación del segundo pico correspondiente al Na en este caso. | 
Figura - 6 |
Las diferencias que se observan al comienzo y al final de ambas curvas son debidas a los efectos de integración y ruido en el registro experimental y por tanto no se considera que deban ser ajustadas en la simulación matemática de los picos.
La concentración de Na encontrada supone que el jarabe analizado tiene unos 23 mg/ml, lo cual es el valor máximo posible admitido en el prospecto del mismo.
AJUSTE DE TRES PICOS
Aunque menos frecuente que el caso de dos picos no resueltos, hay ocasiones en que aparecen tres picos que no están completamente resueltos y que interfieren en la cuantificación de uno o varios de los componentes del pico registrado experimentalmente.
En ese caso la dificultad de descomponer el pico experimental en sus componentes puede incrementarse dependiendo, una vez mas, de la proximidad de los tiempos de elución de los picos implicados entre sí, aunque puede suceder que alguno de los picos esté algo mas separado y aunque no esté resuelto pueda proporcionar información sobre su tiempo de elución, la constante F que relacione los tiempos de bajada y subida o alguna otra que pueda ayudar en el proceso.
Como ejemplo de uno de estos casos se toma un pico que se ha obtenido por electroforesis capilar en placa, con las mismas condiciones que las utilizadas en el ejemplo anteriormente expuesto para el ajuste de dos picos.
Para ello se ha realizado la electroforesis de una muestra de agua del grifo de Madrid, y se han inyectado también en el mismo análisis patrones de K, Na, Li, Ca y Mg , para poder cuantificar estos metales y ver si los valores obtenidos se encuentran dentro de los rangos facilitados por el Organismo Oficial de control de la calidad del agua de Madrid.
Igual que en el caso de dos picos, se traslada el pico que se ha obtenido experimentalmente, expandido 6 veces en la escala de tiempos, a la escala en la que también se construirán los picos simulados matemáticamente.
El tiempo total entre el comienzo y final del tramo transferido es 360s por lo que los picos simulados tendrán que atenerse a dicha escala de tiempos y no a los tiempos indicados en el eje de abscisas de la figura. | 
Figura - 7 |
Se puede ver en la figura - 7 que el pico experimental tiene poca intensidad y presenta dos mínimos uno de los cuales, el primero de ellos, coincide con el tiempo de elución del Ca por lo que podemos suponer que la muestra analizada tiene Ca, entre otros componentes.
En un tiempo anterior, coincidiendo con el tiempo de elución del K aparece un pequeño pico, no representado en la Figura - 7 , bien resuelto y que puede ser cuantificado independientemente con respecto al patrón de K inyectado, como se verá mas adelante.
Se ve en la figura que aunque la base del pico no es completamente plana ni constante, se puede ajustar razonablemente bien a un valor constante que se toma como base de intensidad Cb para la construcción de los picos simulados.
Se construye un primer pico cuyo tiempo de elución se ajuste al del primer mínimo del pico experimental y su tiempo de bajada se adapte a el observado en dicho pico.
También este pico proporciona algo de información sobre la intensidad y el tiempo de subida, aunque esto se podrá quizá modificar posteriormente, cuando se ajusten los otros picos que previsiblemente se encuentran no resueltos en el pico experimental.
| 
Figura - 8 |
El segundo pico que está integrado en la segunda parte del pico registrado experimentalmente tendrá que ajustarse también al tiempo de bajada observado en dicho pico y también tener un tiempo de elución que se pueda ajustar, con ese tiempo de bajada, al observado en el pico experimental.
El valor de intensidad que se asigna a este pico será un valor tentativo que despues podrá corregirse al realizar el ajuste final de los picos simulados.
En la Figura - 9 se ve el resultado de incluir este segundo pico en la simulación, donde se aprecia que la curva de bajada del pico simulado se ajusta bien a la experimental y en la que se comprueba que no es posible conseguir el ajuste completo con solo este pico por lo que se deduce que hay otro pico incluido en el registro, que habrá que añadir a la simulación para conseguir su ajuste completo con el pico experimental.
| 
Figura - 9 |
El tiempo de bajada de este pico se ha elegido ligeramente superior al del primero, lo cual es consecuente con el hecho observado de que los tiempos de bajada dependen de los tiempos de elución, siendo mayores cuanto mayor es dicho tiempo.
De esta forma se consigue también el ajuste de la curva de salida del primer pico, lo que no se puede conseguir si mantenemos un tiempo de bajada igual al anterior.
Al incluir un tercer pico, con un tiempo de elución retrasado 11s con respecto al del segundo y con el mismo tiempo de bajada, se observa que el ajuste ha sido mejorado considerablemente.
La constante F elegida al ajustar el primer pico se ha mantenido igual para los otros dos y parece adecuada, a falta de pequeños cambios en los tiempos de elución y las constantes de intensidad que mejoren el ajuste de las dos curvas. | 
Figura - 10 |
El ajuste conseguido que se ve en la Figura - 10 confirma que dentro del pico registrado experimentalmente hay tres picos, uno de ellos mas evidente que los otros dos que se encuentran solapados dentro de un solo pico aparentemente.
Retrasando a 9s el tiempo de elución del tercer pico con respecto al segundo y retocando ligeramente las constantes de intensidad de los picos y los tiempos de bajada, se observa en la Figura - 11 que el ajuste ha mejorado con respecto al de la anterior figura, aunque como se mencionó anteriormente se comprueba que estas pequeñas variaciones solo producen diferencias en los valores de cuantificación del orden del 10%. | 
Figura - 11 |
Como consecuencia de este análisis se puede concluir que los metales mayoritarios, presentes en la muestra de agua del grifo de Madrid son ; K, Ca, Na y Mg , dado que los tiempos de elución de los dos picos solapados en el segundo pico no resuelto coinciden con los tiempos de Na y Mg , respectivamente y los picos de K y Ca eluyen en sus tiempos correspondientes.
Las concentraciones encontradas para dichos metales son : K - 1.2 mg/l , Ca - 6 mg/l , Na - 5 mg/l , Mg - 1.3 mg/l.
Estos valores están dentro de los márgenes promedio facilitados por la Subdirección de Calidad de las aguas de Madrid en informe de 15 de Enero de 2026 (4), que son : K (0,5-1,4 mg/l) , Ca (6,4-17 mg/l) , Na (5-13 mg/l) , y Mg (1,1-2,4 mg/l).
REFERENCIAS
4.- Subdirección de Calidad de Las Aguas. ANÁLISIS COMPLETO. Canal de Isabel II, 15/1/2026.
RESULTADOS EXPERIMENTALES
Hasta ahora se han tratado ejemplos relativamente sencillos en los que la composición de los picos era conocida previamente y sus concentraciones no presentaban problemas de detección o saturación en el análisis experimental, pero hay algunos casos que pueden complicar algo mas la interpretación de los picos no resueltos y su descomposición en picos individuales que permitan la cuantificación de los analitos de interés.
Tal es el caso de algunas aguas de consumo embotelladas, calificadas de "minerales" por su alta concentración en sales de metales alcalinos y alcalinotérreos, en las que habrá que recurrir a su dilución previa al análisis para evitar la saturación del detector y en las que la disparidad de concentraciones entre unos y otros componentes puede enmascarar los picos de algunos de ellos, de baja concentración, frente a otros mucho mas concentrados.
Los análisis se han realizado por electroforesis capilar en placa con electrolito de fondo de ácido acetilsalicílico, longitud de columna de 10cm, volumen de inyección de 2ul, tensión de polarización de 500v y detección por conductividad eléctrica.
La primera muestra analizada es un agua embotellada, calificada de "mineralización media" en la que se ha hecho una dilución previa de la muestra a 1:4 y cuya curva resultante del análisis se muestra en la figura siguiente.
Se traslada el pico experimental expandido 6 veces en la escala de tiempos, con respecto a la leyenda del eje de abscisas y aunque se ve que no tiene una linea de base plana y constante, se puede ajustar una linea de base de intensidad Cb en el comienzo del pico para construir los picos simulados.
A simple vista se aprecia en la Figura -12 la presencia de lo que parecen dos picos no resueltos, uno de ellos de menor intensidad que el otro.
| 
Figura - 12 |
Se construye un primer pico para ajustar la curva simulada a la experimental, eligiendo el tiempo de elución, el de bajada y la intensidad del pico de forma que se ajuste al primer pico que se observa en el registro experimental.
Se ve en la Figura - 13 que este primer pico simulado ajusta al experimental, y se ha elegido también una constante F que ajuste la subida del pico simulado a la del experimental, lo que se podrá modificar en el proceso posterior de simulación de la curva.
Se ve en la figura que este primer pico tiene poca intensidad y su tiempo de elución coincide con el del Ca, por lo que se puede suponer que hay una cierta concentración de Ca presente en la muestra.
| 
Figura - 13 |
Siguiendo con la construcción de la curva simulada hay que incluir un segundo pico en la simulación, adaptando su tiempo de elución y tiempo de bajada a los observados en la curva experimental
El tiempo de bajada de este pico se ha supuesto igual al del primero y la constante de intensidad y el tiempo de elución se han elegido para adaptarse a la curva experimental.
Se observa que con esas constantes se adapta bien la curva de salida del primer pico y que no es posible simular la curva experimental con solo estos picos por lo que debe de haber un tercer pico incluido en el registro experimental.
| 
Figura - 14 |
Por ello es necesario incluir un tercer pico en la simulación para poder ajustarse al registro experimental, el tiempo de bajada del mismo se puede considerar igual al del pico anterior dada la proximidad de los tiempos de elución, el valor de cuya diferencia habrá que estimar, ajustando la curva simulada al registro.
Al incluir este tercer pico con una diferencia de tiempo de elución de 14s con respecto al segundo y con una constante de intensidad que luego se podrá adaptar mejor, se ve en la Figura - 15 que la curva simulada se va ajustando a la experimental.
Se ve en la figura que se necesita hacer una adaptación posterior de la constante F para que los tiempos de subida de los picos simulado y experimental coincidan.
| 
Figura - 15 |
Cambiando el valor de la constante F y haciendo ligeros cambios en las constantes de intensidad del segundo y tercer pico, se llega a la curva simulada que se puede ver en la figura siguiente.
Al hacer estos cambios, se puede ver en la Figura - 16 que el ajuste entre la curva experimental y la simulada es aceptablemente bueno.
Las diferencias que se observan al comienzo y al final de la curva experimental con respecto a la simulada son debidas a desviaciones de la linea de base en el registro experimental y efectos de ruido en dicho registro, por lo que no hay que considerarlas en la cuantificación de los picos obtenidos.
| 
Figura - 16 |
A la vista de los tiempos de elución obtenidos se concluye que la muestra contiene Ca , Na , Mg y una pequeña concentración de K, menor de 1 mg/l. Los valores de concentración encontrados para dichos componentes son, una vez aplicado el factor de dilución : Ca - 6.4 mg/l , Na - 40 mg/l , Mg - 6 mg/l.
Los valores orientativos de estos componentes que se especifican en la etiqueta de la botella son : Ca - 9 mg/l , Na - 44 mg/l , Mg - 5 mg/l , y no se hace mención a la presencia de K , quizá debido a su pequeña cantidad.
Este caso presentaba una de las posibles dificultades que se pueden encontrar en el análisis de muestras de aguas y se ha solucionado diluyendo la muestra y realizando el análisis sobre la muestra diluida, pero hay otro tipo de muestras en las que el problema es el contrario, tienen muy poca concentración de sales y ya no se puede resolver diluyendo la muestra porque habría que hacer lo contrario, es decir concentrarla, lo cual puede resultar engorroso y ser una fuente de errores.
También se puede recurrir en estos casos a aumentar la ganancia en el registro experimental, pero esto aumentará también el ruido y las posibles derivas del mismo lo cual tampoco resulta deseable para el ajuste de la posterior simulación matemática del registro.
Como ejemplo de una muestra de este tipo, se ha realizado la electroforesis de un agua de consumo embotellada calificada como de "mineralización débil" . Las condiciones del análisis han sido las mismas que las de la anterior muestra, manteniendo el volumen de inyección en 2 ul.
Se traslada el pico experimental obtenido, expandido 6 veces en la escala de tiempos, a la escala sobre la que también se van a construir los picos de la simulación matemática.
Se puede ver en la Figura - 17 que el pico experimental tiene una intensidad muy pequeña, pero al menos tiene una linea de base plana y constante que permite asignar a la misma una base de intensidad Cb para la simulación matemática de los picos.
| 
Figura - 17 |
Se forma un primer pico simulado eligiendo el tiempo de bajada y el tiempo de elución para que la curva simulada se ajuste a la bajada de la curva experimental, como se puede ver en la figura siguiente.
La constante de intensidad del pico se toma de forma que el pico construido se ajuste bien al experimental y posteriormente se podrá adaptar la constante F de relación entre tiempo de bajada y de subida.
En la Figura - 18 se puede ver que la curva experimental no puede simularse solo con este pico, aunque se ajustara mejor la constante F, lo que indica que hay que incluir en la simulación al menos otro pico.
| 
Figura - 18 |
El tiempo de elución elegido para el primer pico coincide con el del Ca, por lo que al tratarse de un agua de consumo es muy probable que los otros componentes incluidos en el pico sean Na y/o Mg, por lo que se incluye en el segundo pico un tiempo de elución con una diferencia de tiempo respecto del Ca consecuente con esta hipótesis.
Con un tiempo de elución separado por 33s respecto del primer pico y manteniendo el mismo tiempo de bajada que el primero se obtiene el pico, al que se asigna una constante de intensidad que posteriormente se puede variar, que se puede ver en la Figura - 19.
Se comprueba también que solo con estos dos picos no se puede conseguir un ajuste aceptable de la curva experimental con la simulada.
| 
Figura - 19 |
Por eso se añade a la simulación un tercer pico, esta vez con una diferencia de tiempo de elución respecto del segundo de 17s , manteniendo igual el tiempo de bajada y con una constante de intensidad que podrá ser posteriormente ajustada.
En la figura - 20 se puede ver el resultado de añadir este tercer pico y todavía resulta dudoso el ajuste de la curva simulada respecto de la experimental.
Se ve que los tiempos de subida de los picos necesitan un ajuste y por ello se actúa sobre la constante F para ver si con ello y algún cambio mas en las constantes de intensidad de los picos segundo y tercero se consigue mejorar el ajuste entre ambas curvas.
| 
Figura - 20 |
Realizando estos cambios se consigue que las dos curvas coincidan aceptablemente bien, como puede verse en la Figura - 21
Se ha conseguido un ajuste aceptable de la curva simulada a la experimental, pero también resulta evidente que en este caso la dificultad es mayor dadas las pequeñas intensidades de los picos que hay que manejar para ello.
Es cierto que pequeñas variaciones en las constantes utilizadas, producen variaciones en las concentraciones calculadas del orden del 10% y al ser concentraciones pequeñas esto no supone variaciones muy grandes en los valores calculados.
| 
Figura - 21 |
Como conclusión del análisis se hace una cuantificación de los picos generados matemáticamente, con respecto a los patrones inyectados, y que corresponden a concentraciones de Ca, Na y Mg presentes en la muestra de agua analizada.
Los valores encontrados son ; Ca - 4.2 mg/l , Na - 1.1 mg/l , Mg - 0.6 mg/l , que son concentraciones muy bajas como ya se anunciaba en la calificación del agua como "mineralización débil" y en acuerdo aceptable con los valores orientativos declarados en el etiquetado de la botella, que son ; Ca - 6.3 mg/l , Na - 1.2 mg/l , Mg - 0.8 mg/l .
Como se esperaba, cuanto menores sean las concentraciones analizadas se estará mas cerca del límite de cuantificación de la técnica de análisis y la exactitud de los resultados será mas cuestionable, si bien en este caso se puede considerar aceptable.
| |
